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Philip Hall 
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Abstract 

We consider the centrifugal instability of the viscous fluid flow between 
concentric circular cylinders in the small gap limit. The amplitude of the 
Taylor vortex is allowed to depend on a slow time variable, a slow axial 
variable, and the polar angle 0. It is shown that the amplitude of the 

vortex cannot in general be described by a single amplitude equation. 
However, if the axial variations are periodic a single amplitude equation can 
be derived. In the absence of any slow axial variations it is shown that a 
Taylor vortex remains stable to wavy vortex perturbations. Furthermore, in 
this situation, stable non-axisymmetric modes can occur but do not bifurcate 
from the Taylor vortex state. The stability of these modes is shown to be 
governed by a modified form of the Eckhaus criterion. 
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In two recent papers labeling [1] and Brand and Cross [2] have 
independently proposed an amplitude equation which governs the slow azimuthal 
and axial evolution of a Taylor vortex in the small gap limit. In this note 
we show that this amplitude equation corresponds to a velocity field which 
necessarily violates the no-slip condition at one of the cylinders. The 
remedy for this difficulty is well-known in hydrodynamic stability theory 
following the work of Davey, et al. [3] and requires the insertion of an 
eigenfunction in the expansion of the disturbance pressure field. The 
presence of this eigenfunction means that the evolution of a Taylor vortex 
cannot be described by a single amplitude equation. 

We shall see that if axial variations are Ignored then it is possible to 
describe the azimuthal evolution of a Taylor vortex by a single amplitude 
equation. However, even this reduced equation differs from the reduced form 
of the equation of labeling, Brandt and Cross. The appropriate amplitude 
equation is discussed in some detail and it is shown that in the small gap 
limit a Taylor vortex is stable to wavy vortex perturbations. Thus the 
evolution equation approach to describe the azimuthal evaluation of a Taylor 
vortex gives results which are not consistent with the classical results of 
Davey, et al. [5] and the available experimental results. The implications of 
this situation will be discussed later. 

We consider then the stability of the flow of a viscous fluid of 
kinematic viscosity v between cylinders of radii R^, + d. The outer 

cylinder is held fixed whilst the inner one rotates with angular velocity 
We define the Reynolds number R and the parameter 6 by 
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where Tq = 3390 Is the critical Taylor number in the small gap limit* The 
velocity components in the radial, azimuthal and axial directions are scaled 
on v/d, and v/d, respectively whilst the pressure is scaled on 

P v^/d^. 

In the small gap limit the basic flow (0, v, 0) has the 

asymptotic form 

v=l-x + 0(6) 
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where x is a radial variable scaled on d. The equations governing the 
stability of this basic flow can be written in the form 


Lu = - + Ql + Tv V + 0(6) 


(2a) 


T J~’^ 9p , J. , 

Lv=-/ = 6 ^+0„ + u-^- + 


T ” 90 ■ ^2 ^ dx 


(2b) 


Lw = - + 0(6) 


(2c) 


9u , y~T f. V 2 9v , 9w 

K— + /-;r 6 ^— +_ = 


9x 


96 9z 


0(6), 


(2d) 
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where z and t have been scaled on d and v/d respectively, 
nonlinear functions Qj^, and are 0(6*^) whilst the operator L 

defined by 


L E 


9 

9t 


+ 6 


V2 /I 

2 





The 

is 

(3) 


Since the Taylor number differs from its critical value by 0(6) we 

Vo 

expect a finite amplitude motion of 0(6 ^ ) and therefore expand u = 

(U,V,W) in the form 


together with a similar expansion for the pressure. We then define the slow 

Vo - 

time scales t — 6 '^t, t — 6t. The details of such an expansion procedure 

V 2 


follow closely those of Tabeling, Brand and Cross and at order 6 
found that 


it is 


_ A(9,T,t)e^®^ 

■iiO 2 


(Uq(x),Vq(x),Wq(x)) + COMPLEX CONJUGATE 
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where a is the critical axial wavenumber whilst (Uq, Vq, Wq) Is the 
velocity eigenfunction corresponding to the critical point on the neutral 
curve. At order 6 it is found that A must satisfy the equation 



M M 

®o ae 3 t 


0 


and labeling has calculated sq numerically and found that sq = .5261. It 
follows from the above equation that A = A($,T) where 



At order 6 the first harmonic and mean flow correction are determined and it 

is at this stage that the difficulty overlooked by labeling. Brand, and Cross 

arises. The mean flow correction at this order is in the 0-direction and we 

denote it by Vj^. However, we see from the equation of continuity that this 

3/2 

mean flow drives a radial mean velocity field of 0(6 ) which we denote 

by uj^. The equation which determines Uj^ is 



and this equation must be integrated to satisfy ® x = 0,1. This 

cannot be achieved unless v^j contains some arbitrary function of 0 and 
t. It is for this reason that the solution given by Tabellng, Brand and Cross 
does not satisfy the no-sllp condition everywhere. The remedy is to allow for 
a pressure eigenfunction in the manner discussed by Davey, et al. [3] and 
DiPrlma and Stuart [6]. Thus the perturbation pressure must be expanded in 


the form 
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= 6^/2 


Pn + ^Pi + 


3/2 


_ 1 /, _ _ 

5 Po(0.t,T) + 


where the relatively large size of the induced mean pressure field is, of 
course, a lubrication effect. The equation for Vj^ now becomes 


a2 

9 V, 


M 




9x 


- - / 2 ^0 1 I .|2 d r,, >, 

90 2 dx ^^0 ’ 


(5a) 


which can be integrated subject to ® x = 0,1. We can then 

substitute for into the equation of continuity to find u^j. The 

condition tht Uj^ should vanish at both x = 0 and x = 1 gives 


- 

J~2 "O 
^0 36^ 


^ 30 % 


(5b) 


where 



Fp(x)dx 


with 


X 1 

F (X) = / U V dx - X / U V dx. 
u Q u u ^ u u 


The equation for is now integrated once and the arbitrary constant which 

appears in the resulting expression for is fixed by insisting that Pq 

be periodic in 6 . The function v^ is then completely determined and we 
find 

■ I IaI^Fq(x) + {^ " x} (6) 


Qq 2tt 2 

where = 3 — / |A| d6. 

0 IT 

Brand and Cross corresponds to 


The amplitude equation found by Tabeling, 
1 9 

v^ = - y Ia| Fq(x) so that the radial mean 
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flow Induced at higher order in their expansions cannot satisfy the no-slip 

3 /2 

condition at both cylinders. At order 6 we find that A satisfies the 


equation 


= T' A + c ^ 
_ A -r c, 2 

3t 3$^ 


- c^ AlAl^ + Cgfl Oq - 3Qq |A|^}A, 


where T' = [t^ - T^] with the order 6 correction to the 
axisymmetric critical Taylor number. The constants Cq, C 3 , C 4 are given by 
Tabeling as 


Cq ~ 26.16, c^ = 2.609, c^ = 40.2 with V 2 ) = 1 


whilst C 5 has been calculated by DiPrinia and Stuart [5] who give 


= -4.76. 


The amplitude equation given by Tabeling, Brandt, and Cross corresponds to 
setting C 3 = 0 in (7). The linearized fom of (7). shows that the non- 
axis 5 anmetrlc mode with wavenumber M is linearly unstable for 


TJ > = C3 


and the finite amplitude axisymmetric mode which bifurcates from T' is 

Ic 


A = A = 
e 


Cl - C3 M ) 


and of course only integer values of M have any physical relevance. The 
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first mode to bifurcate Is the Taylor vortex solution which has M = 0. In 
order to Investigate the stability of (9) we write 


A = Ag + b 


and the linearized equation satisfied by b Is 


= C 3 - (t' - C 3 M^)(b + be^^^^)(l + e) + c. b 

9 1 9$ 


e e 


2tt 


- C3M^] (be^^^+ be^”^)d<|, 


( 10 ) 


where 

3Qq C 5 
e = . 

If we set M = 0 In (10) we can study the stability of a Taylor vortex to ^ 

dependent perturbations. We can then see from (10) that the growth rate of a 

disturbance proportional to icos or Isin M<I> is -C 3 M^ so that, in the 

small gap limit, there is no bifurcation from a Taylor vortex to a wavy vortex 

solution. The non-axisymmetric modes with M ^ 0 are susceptible to the 

Eckhaus-Ben jarain-Feir sideband instability mechanism. Following Stuart and 

DlPrlma [7] it can be shown from (10) that the non-axis ymme trie mode is 

unstable to sidebands with integer wavenumbers yM,( 2 - y)M for 

“1<Y<3,y?'= 1. The non-axisymraetrlc mode with wavenumber M which 

bifurcates from T^ is found to be unstable to such a sideband for 
ic 



< T' < [l + 


(Y+1)(3-y)i^, 
2(14€) J^lc 
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which reduces to Eckhaus^ result TJ < in the limit y 1 with e = 0, 

We see then that the Eckhaus criterion is altered if e ^ 0 so that the 
pressure eigenfunction decreases the unstable regime. In the present problem 
only integer values of yM are physically acceptable so that the non- 
axisymmetric mode is unstable for 


Ic 


< T' < [1 


[4 - 




2(1+£) Ic 


It is of course possible that stable wavy vortex solutions to (7) exist 

but do not bifurcate from Taylor vortex flow. In order to investigate such a 

possibility (7) was integrated numerically using a fully implicit finite 

difference scheme. Several runs were made using difference initial 

distributions for A at different supercritical values of T'. At 

sufficiently large values of t the numerical solutions approached one of the 

equilibrium solutions (9) for some value of M, we found no steady state 

solution other than those given by (9). 

In previous descriptions of non-axi symmetric motion in the small gap 

limit it has become customary to take 2R 6 = T and ^ ~ 6 ^ 

equations of motion. This procedure was followed by labeling, Brand, and 

Cross who then allowed T to differ from its critical axisymmetric value by 

0(e) and considered azimuthal variations with ~ » 1. The 

procedure outlined in this paper has ~ 0(1) but, by taking e = 6 in the 

work of Tabeling, Brand and Cross and inserting the required pressure 

eigenfunction, we can recover (7). Alternatively by taking the limit T' -»■ “ 

3 Vo Vo 

with ~ ^1 * A ~ T' in (7) we recover the corrected form of the 

equation of Tabeling, Brandt, and Cross. Thus the expansion procedures are 

equivalent and it is not necessary in the small gap limit to ignore the 

0 

available experimental results and take » 1. 
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We shall now derive a generalized form of (7) which takes account of slow 
variations of the vortex amplitude in both the axial and azimuthal 
directions. Such an equation has been given by labeling, Brand, and Cross but 
the velocity field associated with that equation does not satisfy the no-slip 

g 

condition at one of the cylinders. We again assume that 0(1) and now 
choose to look for axial variations on the same length scale as those in the 
azimuthal direction. We therefore define 

? = 6 2 


and retain the expansion (4). However we now write 

^ = A(e , .g . .T. ;0 [u^(x), Vq(x), Wq(x)]6^®^ + COMPLEX CONJUGATE, 

and note that in the absence of any pressure eigenfunction the radial mean 
3 /2 

flow at order 6 cannot satisfy the no-slip condition at both cylinders. 

We therefore expand the pressure in the form 

p = 6 + ... + 5 p^(e ,T ,t) + pj^(9 ,C ,T,t) 


+ 5 ^^^2 p2(0,?,T,t) + ••• , 


( 11 ) 


where ~p^y etc. are pressure eigenfunctions introduced in order that 

the radial mean flow satisfies the required boundary conditions. However, it 
follows from (2c) that Pq will drive an axial mean flow of order 6^ if 
8Po 9Pn 3Pl 

•T — * 0. Thus we set -k — = 0 and -s— - = 0 in order that the azimuthal and 

9? oi, dg 


* 0 . 


9? 
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axial mean flows Induced by the disturbance should be comparable. If we 
denote the mean part of Wj^ by Wj^j it follows from (2c) that 


W~ ’ 


so that 


1 ^P2 

2 8? 


(x - x), 


whilst vj^, the order 6 azimuthal mean flow again satisfies (5a) so that 


v„ - - I lAl' - 


9 ^Pn ^ 

2 ^0 rx - X>| 

Tq 80 2 ^ 


The radial mean flow at order 6 Is now driven by both Vj^ and Wj^ and 

satisfies the no-slip condition at x = 0,1 if 


^ + 6i- lAl^ 0 A 


^0 2 


where Qg is as defined earlier. We see that if = 0 then the pressure 

eigenfunctions pg and both vanish whilst if = 0 we recover (5b). 

It is of Interest to note that in the corresponding pressure equation in 

parallel flow stability theory Pq “ ^ 2 * 

The axial mean flow of order £ interacts with the fundamental terms of 
. 1/2 

order 6 to reproduce the fundamental so that the solvability condition at 
3/2 

order o will depend on . After some analysis we find that the 

appropriate solvability condition is 
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3l 


2 2 

-r- A ^ 9 A ^ 3^ A 


+ ic 


3^ A 
7 35 36 


- c. 


AlAl 


9P0 9P2 

“ *^5 30 ^ “ ^*^8 35 


(15) 


where C 3 , c^, and are as defined previously whilst labeling has 

calculated Cg, cy and gives 


Cg = .9837, c^ - .3948. 

The constant Cg satisfies 

a I [j V+(x 3 - x)V„ - i uJCx^ - x)(D5 - a\) + 0„ ujdx 

Cg = J , (16) 

- a^ U„) - vj V„]dx 

where (^q> ^q) function pair adjoint to (Uq, Vq) . We have not 

calculated Cg which is purely real but there is no reason to suppose that it 
is zero. The equations (14) and (15) are coupled so that in general it is not 
possible to describe the axial and azimuthal evolution of a Taylor vortex by a 
single amplitude equation. 

In order to demonstrate how (14) and (15) can be solved let us suppose 
that all the disturbance quantities are periodic in C with wavelength 2L. 
We can integrate (14) once to give 

^= 6 Qo/^|e / Ia| 2 d5' - (5 + L) + F(0,r), (17) 


where F is an unknown function of 0 and t. However, if P 2 is periodic 
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in ^ 
so that 


), (19) 

and if we integrate (19) once more with respect to C and impose the 
condition that periodic in ? we obtain 

F( 6 ,t) = / /{|A|-^ / ur dr}dr d? . ( 20 ) 

“L ™’L ”Xj 

The solution for p^ obtained by this procedure is unique only up to an 
arbitrary function of 0 and t. This nonuniqueness of p« can be seen 

3p2 9 P 2 

from (14) and (15) directly, since only • 5 —— and — ^ — appear in these 
equations . 

We now integrate (18) to give 

8 ? 3Q /f7 L 2 tt 

30^= "IT ’'4 ^ {lAl -^/ IA|2 d0}d?, (21) 

”L 0 

where the unknovm functions of t which results from integrating (18) has 

9P2 9?o 

been fixed so that p^ is periodic in 0 • The functions and can 

then be substituted into (15) to give an equation depending only on A which 
must then be integrated numericallly subject to the periodicity of A in 0 
and C • In the more general case where the perturbation is not periodic in 


- / 
afl J 


Ur dr}d?' + F(0,t 



C such a simplification of (14) and (15) is not possible. 
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CONCLUSION 

We have given a self-consistent derivation of the equations governing the 
evolution of a Taylor vortex flow in the small gap limit. The equations which 
we have derived differ from those derived previously for the same problem. 
Our results for the situation where there is no slow axial variation suggest 
that, unlike the finite gap case, there is no bifurcation to a wavy vortex 
flow from a Taylor vortex flow. Furthermore we have found that non- 

axisymmetric modes are possible stable states in the small gap limit. This 

latter result is again different from the classical results of Davey et al. 

[5] for the finite (but small) gap problem. The expansion procedure used in 
the latter paper allowed the axisymmetric and non-axlsymmetric eigenvalues to 

be split apart by taking yg- » 1 but results were obtained by subsequently 

9 9 

taking 0(1)- The expansion procedure used here has yg" 0(1) 

throughout and at first order in our expansions we have a multiple eigenvalue 

so that a self-consistent expansion can be developed. It is possible that the 

procedure we have used is valid only for very small values of 6 whilst that 

of Davey et al. [5], though not a formally self-consistent gives good 

agreement with experiment at moderately small values of 6. 
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